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Abstract 



The magnetic relaxation of ferromagnetic powders has been studied for many 
years, largely due to its importance to recording technologies. However, only 
recently have experiments been performed that resolve the magnetic state 
of individual sub-micron particles. Motivated by these experimental devel- 
opments, we use droplet theory and Monte Carlo simulations to study the 
time and field dependence of some quantities that can be observed by mag- 
netic force microscopy. Particular emphasis is placed on the effects of finite 
particle size. The qualitative agreement between experiments on switching 
and our simulations in individual single-domain ferromagnets suggests that 
the switching mechanism in such particles may involve local nucleation and 
subsequent growth of droplets of the stable phase. 
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I. INTRODUCTION 



Magnetic recording technologies are among the most important apphcations of perma- 
nent magnets today. A typical storage medium consists of fine magnetic particles suspended 
in an organic binder adhering to a polymer substrate. During the recording process, differ- 
ent regions of the medium are briefly exposed to strong magnetic fields, so that each grain 
is magnetized in the desired direction [Q. Each grain could thus in principle store one bit 
of data, so greater storage density could ideally be achieved by a medium containing many 
small grains than by one containing a few large grains. However, in order to serve as reli- 
able storage devices, the grains must be capable of retaining their magnetizations for long 
periods of time in weaker, arbitrarily oriented ambient magnetic fields. Since experiments 
show the existence of a particle size at which the magnetizations are most stable (see, e.g., 
Ref. 0), there is a tradeoff between high storage capacity and long-term data integrity 
which must give rise to an optimum choice of grain size for any given material. During both 
recording and storage, the relationship between the magnetic field, the size of the grain, and 
the lifetime of the magnetization opposed to the applied magnetic field is therefore of great 
technological interest. 

Fine ferromagnetic grains have been studied for many years, but until recently such 
particles could be studied experimentally only in powders (see, e.g., Ref. 0). This made 
it difficult to differentiate the statistical properties of single-grain switching from effects 
resulting from distributions in particle sizes, compositions, and local environments, or from 
interactions between grains. Magnetic force microscopy (MFM) (see, e.g., Refs. and 
Lorentz microscopy (see, e.g., Ref. |^) now provide means for overcoming the difficulties in 
resolving the magnetic properties of individual single-domain particles. 

The standard theory of magnetization reversal is a mean-field treatment due to Neel 
and Brown [p!0| , p!l| . In order to avoid an energy barrier due to exchange interactions between 



atomic moments with unlike orientations, Neel-Brown theory assumes uniform rotation of all 
the atomic moments in the system. The remaining barrier is caused by magnetic anisotropy 
T^ , either intrinsic or shape-induced. Anisotropy makes it energetically favorable for each 



atomic moment to be aligned along one or more "easy" axes. Buckling, fanning, and curling 
are, like uniform rotation, theoretical relaxation processes with few degrees of freedom and 
global dynamics ]I|JT^ . In all of the above approaches, the lifetime r of the magnetization 
increases exponentially with the particle volume, L'^: 

T oc expiPL^Af) , (1) 

where d is the spatial dimension, i3~^ = ksT is the temperature in units of energy, and A/ 
is the height of the free-energy density barrier to be crossed in the relaxation process. 

However, for highly anisotropic materials there exists an alternative mode of relaxation 
with a typically much shorter lifetime. Small regions of the phase in which the magnetization 
is parallel to the applied magnetic field (the "stable" phase) are continuously created and 
destroyed by thermal fluctuations within the phase in which the magnetization is antiparallel 
to the field (the "metastable" phase). As long as such a region (henceforth referred to as a 
"droplet") is sufficiently small, the short-ranged exchange interaction with the surrounding 
metastable phase imposes a net free-energy penalty, and the droplet will, with high prob- 
ability, shrink and vanish. Should the droplet become larger than a critical size, however. 
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this penalty will be less than the benefit obtained from orienting parallel to the magnetic 
field, and the droplet will with a high probability grow further, eventually consuming the 
grain. The nature of the metastable decay thus depends on the relative sizes of the grain, the 
critical droplet, the average distance between droplets, and the lattice constant, as discussed 
in detail, e.g., in Refs. [p!^-p!7|. 



At this point we wish to emphasize the difference between a droplet and a domain. 
Whereas they are both spatially continuous regions of uniform magnetization, a domain is 
an equilibrium feature whereas a droplet is a strictly non-equilibrium entity. A droplet may 
either grow or shrink, but it will not remain constant. All of the particles referred to in 
this paper are single- domain, which means that in equilibrium their magnetizations will be 
uniform, but as we shall see they may nevertheless decay through the nucleation and growth 
of droplets. 

Because of its simplicity, the kinetic nearest-neighbor Ising model has been extensively 
studied as a prototype for metastable dynamics (see Ref. [0, and references cited therein). 
In particular, square- and cubic-lattice Ising systems with periodic boundary conditions 
have been used to study grain-size effects in ferroelectric switching [|1^,|1^]. A related one- 
dimensional model has been used to study magnetization reversal in elongated ferromagnetic 
particles |^D|-Q. In this article we use the kinetic Ising model to estimate the switching 
field, -ffsw; at which magnetization reversal is thermally induced on experimental timescales 
for given temperatures and system sizes. (Figure |I] is a sketch of switching-field curves 
and the four regions of metastable decay for a typical metastable system with short-range 
interactions.) We also find explicit functional forms for what is actually observed in MFM 
experiments: the time- and field-dependent probability that a grain of given size remains 
unswitched. 

The organization of the remainder of this paper is as follows. In Sec. |T| we describe the 
kinetic Ising model and the Monte Carlo methods we use to study metastable decay. A 
brief summary of the relevant results of droplet theory is given in Sec. |T|. We present our 
simulational estimates of quantities observable by magnetic force microscopy in Sec. In 
Sec. we conclude with a qualitative comparison with some recent experiments and discuss 
some of the implications of this work. 



II. MODEL AND NUMERICAL METHODS 

The model is defined by the Hamiltonian 

n = -Jj2 s^SJ - HL'^m + DL'^m^ , (2) 

where Si = ±l is the z-component of the magnetization of the atom (spin) at site i, J>0 is 
the ferromagnetic exchange interaction, H is the applied magnetic field times the single-spin 
magnetic moment. The sum runs over all nearest-neighbor pairs on a square (generally 
hypercubic) lattice of side L. In this work we do not study the effects of grain boundaries, 
so periodic boundary conditions are imposed. The dimensionless system magnetization is 
given by 

m = L-''Y.Si, (3) 
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where the sum is over all L'^ sites. The lattice constant is set to unity. 

The last term in Eq. @j allows for an approximate treatment of the magnetostatic energy. 
While a more detailed treatment using magnetic dipole-dipole interactions would be more 
physical, this approximation has a significant advantage in computational efficiency. For 
systems with periodic boundary conditions, the first and third terms of Eq. (0) are equal 
when the system size is given by 

(4) 

where (Joo(T) is the surface tension along a primitive lattice vector in the limit L oo. For 
the two-dimensional Ising model, aoo{T) is known exactly The length scale on which 



we would expect a transition from single-domain to multi-domain behavior is approximately 
L/). In this work we set 1^ = 0, which reduces Eq. (0) to the standard Ising Hamiltonian and 
guarantees that the system will be single-domain at equilibrium for all system sizes. Studies 
of systems with D > will be reported elsewhere . 



The selection of the Ising model is equivalent to requiring that the anisotropy constant be 
much greater than any other microscopic energy scale in the system. Since a large anisotropy 
is in fact required for magnetic materials used in magnetic recording media this ideal- 
ization may not be too unreasonable. Simplicity is our main reason for choosing periodic 
boundary conditions and a two-dimensional system, particularly since many equilibrium 
properties of the two-dimensional Ising model in zero field are known exactly |Q and since 



the kinetics of metastable decay has been extensively studied for this model. As a result, our 
model systems more closely resemble ultrathin magnetic films than magnetic grains. A more 
realistic simulation of three-dimensional grains is planned for later study, but we emphasize 
that we expect that droplet theory applies to almost any spin model with high anisotropy. 
Accordingly, equations are written in forms appropriate for arbitrary dimensionality d, even 
though simulations are only carried out for d = 2. 

The relaxation kinetics is simulated by the single-spin-flip Metropolis dynamic with 
updates at randomly chosen sites. This dynamic is realized by the original Metropolis 
algorithm |]2B| and the n-fold way algorithm [^. (For a discussion on the equivalence 
of the dynamics of these algorithms, see Ref. p8[.) The acceptance probability in the 
Metropolis algorithm for a proposed flip of the spin at site a from Sa to — is defined 
as W{sa — Sq) = min[l, exp(— /5Ai?a)], where Ai^^ is the energy change due to the flip. 
The n-fold way algorithm is similar, but involves the tabulation of energy classes. First an 
energy class is chosen randomly with the appropriately weighted probability. A single site 
is then chosen from within that class with uniform probability and flipped with probability 
one. The number of Metropolis algorithm steps which would be required to achieve this 
change is chosen from a geometric probability distribution [^, and the time, measured in 



Monte Carlo steps per spin (MCSS), is incremented accordingly. The n-fold way algorithm 
is more efficient than the Metropolis algorithm at low temperatures, where the Metropolis 
algorithm requires many attempts before a change is made. 

We study the relaxation of the dimensionless system magnetization starting from an 
initial state magnetized opposite to the applied field (m(t = 0) = +1, if < 0). This approach 
has often been used in previous studies, e.g. in Refs. [0,^]. For the temperatures employed 
in this study, the equilibrium spontaneous magnetizations in zero field are close to unity. 
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with 0.95<ms<l. Since the apphed field is negative (and generally small), the equilibrium 
magnetization is assumed to be approximately — mg and the metastable magnetization is 
assumed to be given by mj^s^+rris. We use as an operational definition of the lifetime r of 
the metastable phase the mean first-passage time to a cutoff magnetization m = 0: 

r=(t(m = 0)). (5) 

It has been observed |TB[ that the qualitative results discussed below are not sensitive to 



the cutoff magnetization as long as it is sufficiently less than mg. Our choice of m = as the 
cutoff facilitates comparison with MFM experiments, which are only capable of measuring 
the sign of the particle magnetization. Note that the mean first-passage time is generally 
not the same as the median first-passage time. 



III. DROPLET THEORY 

In this section we briefiy review some elements of droplet theory. We concentrate specif- 
ically on the form of the lifetime within the four regions (shown in Fig. |l|) of metastable 
decay for a typical metastable system with short-ranged interactions. As we will see below 
[Eqs. (H), (H), ([T^)], the exponential dependence of the metastable lifetime of a system with 
short-ranged interactions does not depend exponentially on the system volume, as is the 
case for a system with mean-field interactions [Eq. (|l|)]. For a more complete review, see 



Ref. \17 



At zero applied field, the system enjoys true coexistence between two degenerate equi- 
librium phases with magnetizations m ^ img ~ ±1. This leads to the identification of a 



"Coexistence" (CE) region ||T5|,|T6[ within which the applied field is so weak that metastable 
decay proceeds only through a single, system-spanning subcritical droplet. The field limit- 
ing the CE region is called the "Thermodynamic Spinodal" (i^ThSp) HT^H- '^^^ radius of a 
critical droplet, given by [0 

provides a criterion for estimating the crossover to decay by a single droplet. We estimate 
l-f^ThSpl by setting 2Rc = (L — 1), which yields 

" (L - • 

For a short-range-force system in the CE region the free-energy barrier that must be 
crossed comes mostly from the creation of an interface at intermediate system magnetiza- 
tions. The resulting form for the lifetime for systems with periodic boundary conditions is 

mm 

TiL,H,T)^A{T)exp{2PaLiT)L''-'[l + 0iHL)]} , (8) 

where A{T) is a non- universal prefactor and criiT) is the L-dependent surface tension along a 
primitive lattice vector. Any power-law prefactors in Eq. (|]) are absorbed into this definition 
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of (Jl{T) as 0[(lnL)/L] terms. The 0{HL) term comes from bulk contributions to the free 
energy. 

The region of fields just stronger than i^xhSp is called the "Single-Droplet" (SD) region 
|l5| , p!6[| . In this region the first critical droplet to nucleate almost always grows to fill the 
system before any other droplet has a chance to nucleate. The lifetime in the SD region is 
given in terms of the nucleation rate per unit volume / by 



t{L,H,T) ^ L'^I{T,H) 



(9) 



\H 



i-d 



The nucleation rate /(T, H) is given by 

IiT,H)^BiT)\Hfexp{- 

where B{T) is a non-universal prefactor, and K is believed to be 3 for the two-dimensional 
Ising model and —1/3 for the three-dimensional Ising model Here Ho(T) is given by 



EoiT) + Er{T)H']} , 



(10) 



m 



d 



2m^ 



(11) 



where is a shape- and dimension-dependent constant such that the volume of a droplet is 
given in terms of its radius by V = QR'^. For the two-dimensional Ising model, Q is exactly 



known |T^. We determine Si from a numerical fit. 

In both the CE and SD regions, switching is abrupt, with a negligible amount of time 
being spent in configurations with magnetizations other than img. This phenomenon, in 
which the entire system behaves as though it were a single magnetic moment, is known as 
superparamagnetism |T^J3^. Switching is a Poisson process in both the CE and SD regions, 
so the probability P{m>0) that the magnetization is greater than zero decays exponentially 
with time, and the median switching time is larger than the mean switching time by a factor 
of In 2. The standard deviation of the switching time for an individual grain is approximately 
equal to the mean switching time, r. Because of the random nature of switching in the CE 
and SD regions, their union has been called the "Stochastic Region" [T5|JT6[1 . 

Still stronger fields, by contrast, lead to decay through many weakly interacting droplets 
in the manner described by Kolmogorov [Q, Johnson and Mehl and Avrami p6|-|38|. 
Such decay is "deterministic" in the sense that the standard deviation of the switching time 
is much less than its mean (see the Appendix). The crossover between the SD region and the 



multi-droplet (MD) region has been called the "Dynamic Spinodal" (DSp) ||T5|,|T6[. Since the 
standard deviation of the lifetime is equal to its mean in the Stochastic region, we estimate 
this crossover by the field H1/2 at which 

72 = L 
2 ■ 



v/(t^(m = 0))-r 
For asymptotically large L, Hdsp tends to zero as [|T6| 



(12) 



DSP 



So(T) 



{d+l)\nL 



i/{d-i) 



1 + 



K-1 

rf2 - 1 



+ 



ln(^ 

B 



K 

T- 1 



In 



—0 
d+1 



In(lnL) 

InL 
ln(2fi)" 
d 



1 



(rf+ l)lnL 



(13) 
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however, extremely large system sizes may be required before this scaling form is observed. 
Within the MD region the lifetime is given by |17 



t{L,H,T) 



+ l)ln2 



IiT,H) 



(14) 



where I{T,H) is given by Eq. (|T0|). Here v is the (nonuniversal) temperature-dependent 
radial growth velocity of a droplet, which under an AUen-Cahn approximation p9|-|^ is 
proportional to the applied field in the limit of large droplets: 



iy\H\ 



(15) 



At a sufficiently high field, nucleation becomes much faster than growth and the droplet 
picture breaks down. The crossover to this "Strong-Field" region (SF) has been called the 
"Mean-Field Spinodal" (MFSp) [|l^,|l^]. A conservative estimate for this crossover field is 
obtained by setting 2Rc = l: 



MFSp I 



{d-l)a^{T) 



(16) 



IV. MONTE CARLO SIMULATIONS OF MFM OBSERVABLES 

The droplet-theory concepts of Sec. |ITl| provide a framework for the calculation of a 
number of observable quantities. In this section we discuss calculations of MFM observables 
from Monte Carlo simulations of the two-dimensional kinetic Ising model as discussed in 
Sec. ^ 

The mean first-passage time r to m = can be found at a fixed magnetic field directly 
from Monte Carlo simulations. Finding the magnetic field corresponding to a given lifetime 
requires an indirect approach. From four fields chosen such that t{L,H,T) ^ t a.t fixed 
system size L and temperature T, we use weighted linear regression to estimate the switching 
field Hs^m t, T) such that r(L, i^gw, T) = t. The results, which are shown in Fig. |^ for three 
different temperatures and a range of system sizes and waiting times, are qualitatively similar 
to experimental measurements of the same quantity made on powders (see, e.g., Ref. [Q) 
and on single-domain particles 0. As Fig. |^ illustrates, the maximum value of ifs„(L; t, T) 
occurs within the SD region. The asymptotically nonzero value of Hsw{L]t,T) in the MD 
region is not observed in actual experiments, perhaps because defects and impurities lead 
to heterogenous nucleation in large grains. Lastly, a comparison of parts (a), (b), and (c) of 
Fig. shows the striking temperature dependence of Hsvj for a fixed waiting time. 

From MFM experiments one does not measure magnetization, but rather detects the ori- 
entation of the magnetic poles, i.e., the sign of m. An experimentally measurable quantity 
that can be used to determine the decay mode is therefore the probability that the magne- 
tization is greater than zero, P{m > 0), as a function of field at fixed waiting time or as a 
function of waiting time at fixed field. The particular functional form of P which is predicted 
by droplet theory depends upon the region (CE, SD, MD, or SF) in which P(m>0)~l/2. 

In the CE region, switching is a Poisson process and the decay of {m{t)) has an expo- 
nential form. However, recrossing (returning to an m > phase from an m < phase) is 
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important because the energy difference between the "stable" and "metastable" phases is 
so smalL This is most easily seen at H = 0, where the two phases coexist and 



P(m>0) = [exp(-2t/r) + l]/2 . 



(17) 



The factor of two in the exponential comes from the fact that r is the mean first-passage 
time to m = 0, not to m = — mg. Figure ^(a) shows MC estimates of P(m>0) as a function 
of time, and Fig. ^(a) shows MC estimates of P{m > 0) as a function of field in the CE 
region. 

Because the energy difference between the stable and metastable phases is large in the 
SD region, the probability of return to m>0 from m<0 is negligible. As a result, P(m>0) 
decays exponentially in time, and we specifically obtain 



P(m>0) = exp 



-t 



t{L,H,T) 



(18) 



where t{L,H,T) is given by Eqs. (D and (|T0|). Note that P(m > 0) ^ ^ ^ = 1/e ^ 1/2 
because the mean of an exponential distribution is not the same as its median. In Fig. |^(b), 
Monte Carlo estimates of P(m>0) as a function of t at fixed H, T, and L are seen to agree 
well with an exponential time decay. Figure §(b) shows Monte Carlo estimates of P as a 
function of H at fixed t, T, and L together with fits based on Eqs. (^), (0), and ([T^). The 
fit shown in Fig. ^(b) comes from a least-squares fit used to determine B{T) and Hi(r) in 
Eq. (|10|) H. The fit is seen to agree excellently with the MC data. 

The form of P(m > 0) in the MD region is somewhat more complicated and is derived 
in the Appendix. Both the switching field and the lifetime have asymptotically Gaussian 
probability distributions. Expanding in time around t = r we find 



P(m > 0) 



2 



1 + erf 



t-T 



where r is given by Eq. (0) with H = Hsw Here At is given by 

^ _ (Hgsw|)-^/V(-^+^)/^ 

and ipd is given in Table |I[ Likewise, an expansion in field around ifgw yields 

1 



(19) 



P(m > 0) 



where 



Ah = (rz/|ifs- 



\d/2 



d + K 



-1 



\-d 



'l-d)Eo + {3-d)E,Hl 



(20) 



(21) 



(22) 



Here i^sw is the switching field corresponding to the time r. 

Valuable dynamical information can be found by fitting the MC data to Eqs. (|T^ and 
(pip. From the fit to Eq. (|l^) shown in Fig. |^(c) we can find v and r at fixed H directly 
(Tab. |lD. In order to find v from Eq. ( ^Tj) it is first necessary to find Si. We therefore fit 
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43 1 . The resulting estimates for 



MC estimates of r to Eqs. (|T0|) and (|T4D to find B and Si 
u are consistent with the estimate from Eq. ( |T9|) (Tab. H). Note that for the range of fields 
shown in Fig. ^(c) the radial growth velocity is less than one lattice constant per MCSS, 
which is to be expected from a single-spin-flip dynamic such as we are using. 

Knowledge of B, Hi, and u also allows us to find an alternative estimate for the Dynamic 
Spinodal. The crossover between the SD and MD regions can be expected to occur where 
the mean waiting time for a droplet to nucleate is approximately equal to the time required 
for a nucleated droplet to grow to volume L'^/2. This leads us to 



L'^/(r, H)] = iu\H\)-^ \{2n)-^/'^L - R, 



(23) 



where Rc is the critical droplet radius given in Eq. (|j). As shown by the dot-dot-dashed 
curve in Fig. H(c), this estimate for the Dynamic Spinodal is qualitatively similar to the 
estimate given by Eq. (|T^ but is quantitatively different. This is not surprising, since the 
Dynamic Spinodal is a crossover region, and its precise location depends on the manner in 
which it is defined. 

It is also possible to approximately count the number of critical droplets which nucleate 
during the lifetime of the metastable state. By multiplying the nucleation rate by the volume 
of the system and the lifetime [from Eq. (H)] we find 



Narop = L''l{T,H) 



(rf+ l)ln2 



I{T,H) 



1 

' d+l 



(24) 



Figure ^ shows a typical realization of magnetization reversal in the MD region. From 
Eq. (I), at T = 0.8Tc and if = —0.2 J a critical droplet consists of approximately 48 overturned 
spins. According to Eq. (p4D, approximately 18 critical droplets should have nucleated in the 
system shown in Fig. |^ by the time m = 0, which is in rough agreement with visual inspection 
of the figure. However, it is difficult to confirm this number with precision directly from the 
"snapshots" since i) the droplets may partially or completely overlap, ii) the droplets may 
have irregular shapes, iii) there are many subcritical droplets, and iv) only a few times are 
displayed. 

As a cautionary note, we analyze data taken in the SF region using equations appropriate 
for the MD region. As Fig. ^d) and Fig. ^(d) show, Eqs. (|1^) and ( pTD can be fitted to the 
MC data for P{m>0) quite well. This is because in the deterministic region, the observed 
magnetization is essentially an average over the magnetizations of regions too distant to 
be correlated, so that by the Central Limit Theorem it tends to a Gaussian distribution as 
L^oo. Interestingly, both the value of i' obtained from the fit in Fig. |^(d) [z/ = 0.875(9)] and 
the value from the fit in Fig. §(d) [z/ = 0.82(2)] are in fairly good agreement with the values 
derived in the MD region (Tab. However, the resulting radial growth velocities (f ~ 3 
lattice constants per MCSS) are too large to be meaningful in a single-spin-fiip dynamic. 
Furthermore, whereas Eq. ( [T^ ) gives a good fit for r in the MD region, it fails when applied 
to the entire deterministic region. Since it will be extremely difficult to determine whether 
the radial growth velocity is reasonable in experimental systems (due to uncertainty about 
the appropriate classical dynamic and attempt frequency), droplet theory should be applied 
only if both Eq. (|I4D correctly gives the field-dependence of the lifetime and H <Hmfsp- 
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V. DISCUSSION 



Due to the importance of magnetic recording technologies in modern society, magnetic 
relaxation has been a subject of study for many years. An early mean-field theory for this 



process was presented by Neel and Brown ||10| , pn]| almost fifty years ago, and the most 
popular current models remain mean-field [p],p!3|. This is understandable, since mean-field 
models are much simpler to solve than models with local dynamics that describe spatial fluc- 
tuations. However, a great deal of research has been conducted on non-mean-field descrip- 
tions of met ast ability, particularly for the simple prototype model of anisotropic magnetism, 
the Ising model. These studies make it possible to obtain a more detailed understanding of 
magnetic relaxation. 

In this paper we have discussed the magnetic relaxation of single-domain magnetic par- 
ticles in terms of the droplet theory of metastable decay. We have performed Monte Carlo 
simulations of the metastable decay of two-dimensional Ising systems, and shown that these 
simulations yield switching fields which vary with system size in a manner which is quali- 
tatively similar to experiments conducted on real single-domain particles (Fig. |^). We have 
also used droplet theory to determine the analytical forms of the probability that the mag- 
netization remains unswitched, and fitted these forms to our simulation data (Figs. ^ and 

The methods used in carrying out these fits are appropriate for use with experimental 
data obtained from magnetic force microscopy. 

Because of their two-dimensional nature, the systems we have simulated with MC ac- 
tually serve better as models of ultrathin islands of magnetic material on a nonmagnetic 
substrate. (For a review of experimental observations of the magnetic states of ultrathin 
films, see Ref. [0.) The symmetry of our lattice and the use of Ising spins require the 
modeled material to have strong perpendicular anisotropy. Several such ultrathin films 
have been observed to have critical exponents consistent with the universality class of the 
two-dimensional Ising model (see e.g. Refs. [p5|-^). It would be interesting to compare mi- 
croscopic observations of magnetization reversal in these films with the kinetic Ising model 
described in this paper. 

Because films have macroscopic lateral extents (L^ 1), it is difficult to study any kind of 
finite-size effects in them. One way to impose a finite size to a film is to grow it epitaxially on 
a stepped surface vicinal to a simple facet so that the width of the step is of a microscopically 
moderate size (see e.g. Ref. |Q and references therein). Of particular interest is Fe(llO) 
on W(llO), since at submonolayer coverages the iron can be grown on step edges of the 
tungsten surface and since the resulting film belongs to the two-dimensional Ising model 
universality class . For such semi- infinite systems Monte Carlo simulations such as have 
been carried out here are not possible. Sophisticated techniques in statistical mechanics 
involving transfer matrices have, however, been successfully applied to the two-dimensional 
Ising model in strip geometry pT| , p2[ and to a quasi-one- dimensional Ising model In 



both cases, it was found that concepts from droplet theory are useful in describing the 
metastable decay. 

It is important to understand the relaxation process in the multi-droplet region clearly, 
both because it is relevant to ultra-thin films and some single-domain particles and because 
of controversies which have arisen over the mechanism of magnetic relaxation [Q]. First of 
all, we re-emphasize that droplets are distinct from domains in that they are non-equilibrium 
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entities. Second, our model systems do not include defects or impurities to serve as nuclei. 
In fact, by using periodic boundary conditions and the initial condition m(0) = +l, all sites 
are completely equivalent; this symmetry is spontaneously broken only by local thermal 
excitations. Third, our model system is not a spin glass even though the time-dependence 
of the magnetization may occur as a stretched-exponential function [Eqs. ([Al| ) and (|7^)]. 
Finally, the simulations we have carried out are for an individual particle with a well-defined 
orientation. This distinguishes it from some treatments of magnetic viscosity, e.g. Ref. |E3 



This last point also emphasizes the importance of experiments like MFM which resolve the 
magnetic state of isolated individual ferromagnetic particles. 
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APPENDIX: 



Avrami's Law [p4| -|38[| gives the volume fraction of the metastable state (or equivalently, 
the magnetization) for systems in which droplets nucleate with a rate (per unit volume) 
I and grow without interacting, although they may overlap. The time-dependent system 
magnetization {m{t)) is given by |P^-P 



(m(t)) ^ ms [2{4>{t)) - 1] 



where 



(0(t)) = exp 



-In 2 



T 



(Al) 



(A2) 



is the volume fraction of the metastable phase and r is given by Eq. ([14|) . The factor In 2 
is necessary to fulfill the requirement m(r) =0, in agreement with the definition of r as the 
first-passage time to m = [Eq. (|^)]. 

The nucleation of droplets is random in both space and time, so only the mean magne- 
tization is given by Eq. (|A1|). The variance in the magnetization is related to the two-point 



correlation function, which has been calculated by Sekimoto ||56| for circular and spherical 
droplets. Specifically, if the system is described by a two-valued field u{r, t) where 



u{r,t) 



1 if the position r is within the metastable phase at time t, 
otherwise. 



(A3) 



then the connected two-point correlation function is given by 
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(u(x,t)M(x + r,t)) - (u(x,t))2 

_ f (m(0, |exp [/t;'^t'^+^^rf (r/2i;t)] - l} , r < 2vt 



where r = | r | , and 



0, 



cos y — 2yJl — y'^ + y In 



r > 2t;t 



(A4) 



IT 



^,(y) = :-(l-yf(l + y) . 



Since is simply an average over u{r,t), with 

(0(t))=L"'^|Kr,t))dr 
the variance of the magnetization is given by 



Var[m(t)] ^ AmlL-'^ j [(m(0, t)M(r, t)) - (M(0,t))^ 



dr 



(A5a) 
(A5b) 

(A6) 
(A7a) 



r 







{u{0,t)u{r,t)) - {u{0,t))' 



dr 



Amldn{(P{t))'' 



Amldn{(f){t)) 

j 1 rl 
X < , 

d Jo 



/-^{exp -l}di/ 



+ y'^-iexp ^ ^ i^-j ^a{y) dy\ . (A7b) 



The total system magnetization at time t is thus seen to be essentially an average over on 
the order of {L/2vtY independent samples. As a result of the Central Limit Theorem, for 
sufficiently large L, the distribution in magnetization is well approximated by a Gaussian 
distribution. Introducing the standard deviation in the magnetization, 



a{t) = y'Var[m(t)] 



(A8) 



(not to be confused with the surface tension ai which is discussed in the text), it is easy to 
calculate the probability that the magnetization is greater than zero: 



P{m > 0) ^ / 
Jo 



a{t)V2TT 



exp 



2cr^{t) 



d/i 



1 + erf 



ait)V2, 



(A9) 



This probability will differ from zero or one only if {m{t))/(x{t) ^ 0. Making a Taylor 
expansion in time for {m{t))/a{t) we find 



P(m > 0) 



1 

2 



1 + erf 



t - T 



(AlO) 
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where At is given by 



(All) 



where u is given by Eq. (|T5D, r is given by Eq. (0) with H = Hsvj, and tpd is the constant 
given by 



1 /•! 



/-lexp[(rf+l)ln2^,(y)/^]] 



-1/2 



(A12) 



Table | hsts numerical evaluations of ipd- Likewise, a Taylor expansion of {m{t))/a{t) in 
field yields 



P(m > 0) 



where 



Ah = (rz/|i/sw|) 



d/2 



d + K 



t 2 



1 + erf 



a — 

L-'^/^Ah 



(A13) 



-1 



;i - Ho + (3 - d) E,Hl ■ (A14) 



Here ifg^ is the switching field corresponding to the time r. 
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FIGURES 



FIG. 1. The relationship between the applied field H and system width L for various fixed 
lifetimes (solid curves) in a typical metastable magnetic system. Four regions are distinguished 
by differing decay processes: the Coexistence region (CE), the Single-Droplet region (SD), the 
Multi-Droplet region (MD), and the Strong-Field region (SF). The CE and SD regions are separated 
by the thermodynamic spinodal (dotted curve). The SD and MD regions are separated by the 
dynamic spinodal (dash-dotted curve). The SF region is separated from the other regions by the 
mean- field spinodal (dashed curve). 



FIG. 2. The relationship between the applied field H and system width L for various fixed 
lifetimes (solid curves) and temperatures as calculated through kinetic Ising model simulations. 
The thermodynamic spinodal (dotted curve) is calculated from Eq. (0). The mean- field spinodal 
is approximated by (dash-dotted curve), calculated from Eq. (p^). The mean-field spinodal 
(calculated from Eq. ([l6|)) is not shown, since for each of the three temperatures it occurs at a field 
much larger than the displayed range, a) A;bT' = 0.65J. b) k-QT = 1.?)^,]. c) k-QT = {).Sk-QTc^l.Slf)'if>J. 
The estimate for the Dynamic Spinodal given by Eq. (^) is shown by the dot-dot-dashed curve. 



FIG. 3. P{m > 0) vs. the time t in MCSS for the kinetic Ising model at a temperature of 
T = O.STc- Note the different time scales, (a) L = 6 and H = OJ, in the coexistence region, (b) 
L = 20 and H = -0.105J, in the single-droplet region, (c) L = 100 and H = -0.34725J, in the 
multi-droplet region. The solid curve is a fit to Eq. (|l9|). (d) L = 30 and H = —3. 5 J, in the 
strong- field region. 



FIG. 4. P{m > 0) vs. the applied field H for a kinetic Ising system at a temperature of 
T = 0.8Tc. (a) r = 4670 MCSS and L = 6, in the coexistence region. The horizontal line is 
(1 -|- e~^)/2. (b) r = 914 MCSS and L = 10, in the single-droplet region. The solid curve is a 
least-squares fit of Eq. (18) to the MC data. The inset shows the fitted curve over a wider range of 
fields, (c) r = 40.7 MCSS and L = 30, 100, and 300, in the multi-droplet region. The solid curves 
are fits of Eq. (|l|) to the MC data. The dashed curve is the fit of Eq. (|2l|) to the MC data for 
L = 100, with the width of the switching region reduced by a factor of (300/100)-'^/2 compare 
with the MC data for L = 300. (d) L = 30 and r = 0.8 MCSS, in the strong-field region. 



FIG. 5. Spin configurations showing the nucleation and growth of several droplets in a typical 
realization of magnetization switching in the MD region. Here L = 120, H = —0.2J , and T = 0.8Tc. 
The times shown are at a) t = 50 MCSS, b) t = 114 MCSS, and c) t = 160 MCSS. Grey squares are 
"up" spins and black squares are "down" spins. 
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TABLES 



TABLE L Constants used in Eqs. (19) and (21) to calculate the width of the switching region as 
a function of time and as a function of field, respectively. The derivation is given in the Appendix. 



d 




2 


0.5628545 


3 


0.5109029 
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TABLE II. Results of fits of Eqs. (||) and (|2T|) to MC data for P{m > 0) at T = O.STc. 



Estimated quantities are given with their statistical uncertainties. 



L 


Method 


T [MCSS] 


Hsw 1 J 


vJ [1/MCSS] 


30 


P vs. H 


40.7 


0.34684(2) 


0.96(1) 


100 


P vs. H 


40.7 


0.347253(2) 


0.92(2) 


100 


P vs. t 


40.8(4) 


0.34725 


0.97(1) 
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